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OTHER VACUUM SOLUTIONS WITH 



ALIGNED POLARIZATION 



The solutions of Khan and Penrose (1971) and of Szekeres (1970, 1972), 
discussed in Chapters 3 (and Section 8.2) and 9, describe the coUision of 
plane gravitational waves with aligned linear polarization. It is convenient 
to consider separately in this chapter a number of other exact solutions 
that satisfy this same condition. 

10.1 A general method 

In situations in which the approaching waves are linearly polarized, and 
their polarization vectors are aligned, it is possible to put W = globally. 
In this case the line element (6.20), and the field equations (6.22a-/) take a 
particularly simple form. Equation (6.22a) can immediately be integrated 
to give (6.24) and there is only one main equation, namely (9.3), which is 
a linear equation in V. For any given solution of this equation, a function 
M can always be found satisfying the remaining equations (6.225, c,/), 
although such a function may not satisfy the required boundary condi- 
tions. 

The integral (6.24) involves two continuous functions f{u) and g{v) 
that are monotonically decreasing for positive arguments. As suggested 
by Szekeres (1972), it is therefore possible to use these as coordinates 
in region IV, although care has to be taken at the boundaries and it is 
not possible to extend these coordinates into the prior regions II and III. 
Using (7.8), the line element in the interaction region can thus be written 
in the form 



Vf + 9 



dfdg-{f + g){e 



:W + e-^d|/2), 



(10.1) 



and the main equation (9.3) becomes 



2{f + g)Vfg + Vf + V, = 0. 



(10.2) 



which is the well-known Euler-Poisson-Darboux equation. 
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For any particular solution of (10.2), the new function S can be found 
by integrating equations (7.9), which now become 

Sf = -Uf + 9)Vf, S, ^ -lif + g)V^\ (10.3) 

In order for any particular solution to be appropriate to describe collid- 
ing plane waves, the function S must satisfy the conditions described in 
equations (7.10) to (7.13). 

Using this notation, the components of the Weyl tensor are given by 
the expressions 

*4 = (2^// + (TT^^/ - (/ + ^)^/) • 

It can thus be seen that, although (10.2) is linear so that different solutions 
for V can be superposed, the associated gravitational waves cannot be 
simply superposed. In addition, the scale factors A and B satisfying 
(6.13) and (6.14) must contain components of S, and from (10.3) it can 
be seen that these components are also non-linear in V. 

Once a solution of (10.2) and (10.3) describing the interaction region 
IV is obtained, the approaching waves in regions II and III can immedi- 
ately be deduced. For example, in region II, we simply put 9 = \ and 
assume that / takes the form / = 1 — (citti)"^ -|- . . . . Then, removing 
the constant from 5', it is convenient to put 

g-M(u) ^ f ^-S{u) ^iQ_5) 

nicisjl + f 

and the approaching wave is then given by 

*4 = ^niCiVfe^d{u) 

- Inlcle"^ (2(1 + f)Vff + - (| + ffVf) Q{u). (10.6) 



The approaching wave in region III is obtained in exactly the same way 
in terms of the function g{v). 

In order to obtain exact solutions, the first step is to obtain a general 
class of solutions of (10.2). Szekeres (1972) has indicated how to integrate 
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this equation using Riemann's method, and has given a general solution 
expressed as a line integral involving a Legendre function of order — |. 
This will be described in Chapter 14. However, because it is very difficult 
in practice to evaluate the integrals which this method involves, this does 
not turn out to be a convenient method for obtaining explicit solutions. 

As an alternative approach, it may be observed that attempting to 
solve equation (10.2) by separating the variables leads to the solution 

V = A{a-f)-'/'ia + g)-'/' (10.7) 

where A and a are arbitrary constants. Thus a general class of solutions 
can be obtained by considering 



(10.8) 



for arbitrary sequences of constants cr^ and Aj. The particular decom- 
position (10.8), however, does not turn out to be particularly convenient 
in the construction of explicit solutions. We will therefore proceed by 
changing the coordinates. 

Before introducing a new set of coordinates, it may be noted in pass- 
ing that Feinstcin and Ibariez (1989) have considered an alternative co- 
ordinate system and have expanded a general solution in a different way 
involving Bessel and Neumann functions of zero order. Their approach 
will be described in Section 10.7. 

The appropriate technique is to consider transformations of the basic 
functions / and g which are here treated as coordinates. It turns out to 
be convenient to put 

/= icos(V' + A), g= ^cosiip- X) (10.9) 

where and A are considered as time- like and space- like coordinates, 
which may then be rescaled by putting 

t = sm'i/;, z = smX. (10.10) 

With this, it may be noted that 

f + g=Vl-tWl-z\ f-g = -tz (10.11) 
and t and z can be expressed in terms of / and g by 



-f\/'2+9+\/'2-9^l + f 
= \[f^f\fi+9 - ^\^9^\^f- 



(10.12) 
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The boundaries of region IV are now the hypersurfaces on which t = ±z, 
and the focusing hypersurface in this region occurs when t = 1. Thus, 
to correspond to the interaction region, the coordinates must satisfy the 
inequahty \z\ < t < 1. 

Using these coordinates, the main equation (10.2) becomes 

{{l-t')Vt)^^-{{l-z')V,)^^ = 0. (10.13) 

A general class of solutions of (10.13) can now be obtained by con- 
sidering variable separable solutions of the form 

V = T{t)Z{z). (10.14) 

With this, (10.13) reduces to the pair of Legendre equations 

(1 - t^)Tu - 2tTt + n(n + 1)T = 

, ^ ^ ^ (10.15) 

{1- z^)Z^^-2zZ^ + n{n + l)Z = 

and a class of solutions of (10.13) can be expressed as a sum of products 

= ^ {anPn{t)Pn{z)+qnQn{t)Pn{z) + PnPn{t)Qn{z) + bnQnit)Qn{z)) 
n 

(10.16) 

where Pn{x) and Qn{x) are Legendre functions of the first and second 
kinds respectively, and Qn, Pn and bn are series of arbitrary constants. 
In general, of course, Legendre functions of non-integer order may also be 
included. 

The Legendre functions of integer order are well known, but it may 
still be appropriate to note just the first few. 

Po{x) = l Qo(a.) = ilog([^) 

Pi(x) = X X /! + x\ 

Qi(x) = ^log -1 (10.17) 



P,ix) = 1(5.3 - 3.) Q^^^^ ^ 1(3^2 _ ^1±|^ 



2^ 



It may easily be shown that the solutions of Khan and Penrose, and 
of Szekeres are included in this class. The Khan-Penrose solution uses 



V = -2Qo{t)Po{z), 



(10.18) 
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and the Szekeres solution (9.4) generalizes this to 



V = -{ki + k2)Qo{t)Poiz) - {ki - k2)Po{t)Qo{z). (10.19) 

It may be noted that the Legendre functions of the second kind Qnit) 
are all singular when t = ±1. The singularity t=l occurs when /+gf = 0, 
and thus can be seen to be associated with the focusing singularity in 
region IV. In fact, in order to satisfy the boundary conditions, at least 
one Legendre function of the second kind must be included in the solution 
of (10.13). This may easily be demonstrated by substituting only the 
products of Legendre functions of the first kind into the condition (7.15), 
which can not then be satisfied. 

Another solution of (10.13) can be obtained by considering separable 
solutions involving a sum rather than the product (10.14). This leads to 
the solution 

V = -|alog(l - 1^) - ^alog(l - z'^) (10.20) 

where a is an arbitrary constant. This can immediately be seen to be 
the obvious solution V = aU which, on its own, gives rise to the Kasner 
metrics as will be shown in the next section. This term may be added to 
the sum (10.16). 

A further solution of (10.2), or (10.13) is given by 

V = dcosh"^ ( ""^J'^ ] (10.21) 
V f + 9 J 

where c and d are arbitrary constants. In fact a series of solutions of 
this type can be used having different values of c and d. Thus a different 
representation of the solution of (10.13) can be written in the form 

^ = E*cosh-( ^-;f_^ ). (10.22) 



Alternatively terms from the sum (10.22) may also be added to those 
of (10.16). However, as will be shown in (10.68), some of these terms may 
also be expressed in terms of products of Legendre functions, and so are 
already included in (10.16). It may be noted in passing that particular 
terms of this type have been included by Feinstein and Ibahez (1989) in 
the class of solutions that will be described in Section 10.7. 

Having obtained an expression for V as any combination of the terms 
(10.16), (10.20) and (10.22), it is then necessary to integrate (10.3) to 
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obtain S. These equations may conveniently be rewritten in terms of the 
coordinates i/j and A defined by (10.9) in the form 

2sin(^ + A) " 

COS COS \ 2 

- Sx = . T\i^i> ~ ^a) • 

2sm(7/; — A) 

However, V is now so general that no complete integral for S has yet been 
found. Instead, we proceed by discussing the particular cases that have 
been obtained. 



10.2 The non-singular 'solution' of Stoyanov 

We may now consider a paper of Stoyanov (1979) in which he claimed to 
have obtained a solution without singularities. His method was to look 
for a regular solution of the field equations in region IV, and then to 
obtain the global solution simply by requiring that the metric coefficients 
be continuous across the boundaries. 

The solution he presented for region IV has the line element 

ds^ = 2(1 ± u ± w)("'-^)/2dudw - (1 ± u ± vy-^'dx^ -{l±u± v^+^dy^ 

(10.24) 

where a is an arbitrary constant. This clearly uses 

f=^±ue{u), g=^±vQ{v) (10.25) 

which does not satisfy the boundary conditions (7.3) which require that 
/ and g must be smooth functions. In fact, it has been shown by Nutku 
(1981) that the discontinuities in the derivatives of / and g on the bound- 
aries of region IV indicate the presence of an infinite discontinuity in the 
Ricci tensor on these hypersur faces. Thus, although (10.24) is a vacuum 
solution inside region IV, it cannot be a global vacuum solution describing 
the collision of purely gravitational waves. 

It may also be observed that, with the positive signs in (10.24) and 
(10.25), / and g are increasing functions that are inconsistent with (7.13). 
Such a possibility can only arise if the impulsive components of the mat- 
ter tensor occuring on the boundary of region IV have negative energy 
density. It is in this way that the singularity has been removed. It is 
only with the possibility of the presence of matter with negative energy 
density that the focusing effect of colliding waves can be avoided. 
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Nutku (1981) also pointed out that the hne element (10.24) is the 
well known cosmological solution of Kasner (1921). He gave the explicit 
coordinate transformation by which the line element (10.24) becomes 

ds^ = dt^ - t'^P'idx^f - t^P-'{dxY - t^P'idx'^f, (10.26) 

where the coordinate t is not the same as that used elsewhere in this 
chapter, and 

2(1 -a) 2(1 + a) - 1 

which clearly satisfies the necessary conditions 

Pi+P2+P3 = 1, pI+pI+pI = 1- (10.28) 

The Stoyanov solution has been obtained with the solution of (9.3) 
given by 

V = aU 

= -a\ogif + g) (10.29) 
= -|alog(l - t^) - |alog(l - z^) 

which is the solution (10.20). With this, the remaining equations in (6.22) 
may be integrated to give 

M = - log(c/V) + 1(1 - a2) log(/ + g) (10.30) 

which clearly cannot be continuous across the boundaries of region IV if 
the junction conditions (7.3) are satisfied. 

It must be concluded that the above solution cannot be interpreted in 
terms of an interaction between plane gravitational waves. Certainly, it is 
not a counterexample of a solution for colliding waves without singularity. 
However, in region IV, it is the well known Kasner solution which features 
regularly in discussions of solution generating techniques. In fact, it turns 
out that this solution can be used as a 'seed' from which other physically 
acceptable solutions may be derived. It will therefore be referred to again 
in later sections. 

This solution in the case when a = 1 has been further investigated 
by Taub (1988a), looking particularly at the properties of the distribution 
valued curvature tensor. 
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10.3 The solution of Ferrari and Ibanez and Griffiths 

Now consider the case with the solution of (10.11) given by 



V = -2aQoit)Po{z) - 2bQ2{t)P2{z) 



—a log 



i±^)-^(3.^-l)((3t^-l)log 



1+t 
1-t 



6t 



(10.31) 

This solution was first presented by Ferrari and Ibanez (1986, 1987a) in 
the case when a = 1 and b = —2. The solution with general parameters^ 
was published by Griffiths (1987). The Ferrari-Ibanez (1986, 1987a) solu- 
tion can be seen to be a modification of the Khan-Penrose solution, while 
the more general case is a generalization of the Szekeres solution (10.19) 
with ki = k2 = a. 

With the solution (10.31), the remaining equations (10.23) may be 
integrated to give 

S = const + (a + bf log{t^ - z^) - (a + bf log(l - t^) 



3{a + b)b 
962 



{l-z^)tlog 

1 



1 + t 
1-t 



2z' 



{1-zr 



il-t^)il-9t^)loe{\±l 



(10.32) 



1+t 
1-t 



+ 9t^-4: 



-^(l-.2)(l-t2) 



t^og' 



1 + t 
1-t 



- 4t log 



1 + t 
1-t 



+ 4 



It may be noticed that S contains the term 

(a+bflogit^-z') = lia+bflog[16C^-f)il + f)il-g)C^+g)] (10.33) 

which includes the terms (7.10) that, with (7.8), are required to ensure 
that e~-^ is continuous across the boundaries of region IV. If the leading 
terms in the power series for / and g take the form 

/ = I - {clur + ..., g=\- {C2vr + . . . , n > 2, (10.34) 

then the required boundary conditions are satisfied if 

{a + bf =2{l-l/n). (10.35) 



^ This solution with a + h = + 1 has subsequently also been presented by 
Li (1989). 
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Figure 10.1 Wave profiles of approaching waves in the Ferrari-Ibanez-Griffiths 
solution with n = 2. These have been obtained from (10.6), but the scale factor 
contained in the bounded part of S has been removed. Profiles are shown for 
the special cases in which b = 2/3 or - 4/3, when the wave front is continuous. 
The more general situation is represented by the cases in which b takes the 
values 1 and - 2. 

By extending this solution into regions II and III and putting ci = 
C2 = 1, it can be seen from (10.6) that the approaching waves contain an 
impulsive wave component ^(2a — b)d{u) only if n = 2 and a + b = ±1. 
This case includes the Khan-Penrose solution for which 6 = and the 
Ferrari-Ibanez (1987a) solution in which a = 1, b = —2. However, with 
b ^ 0, the impulsive wave is followed by another wave component. Some 
profiles for approaching waves of this type are illustrated in Figure 10.1. 

It can also be shown using (10.6) that the wave front of the approach- 
ing wave in region II behaves as 

*4 = ^(2a-6)«(i/2-i/-)5(„) 

2 (10.36) 
+ ^(2a - b) ((2a - bf - 4) u^'/^-^Qiu) + 0(«2). 

Thus, unless 6 = 2a or 2a — 2, the wave front may have a distributional 
amplitude. It is unbounded if 2 < n < 4 and has a step if n = 4. The 
wave front is smooth if n > 6. Scaled profiles for the approaching waves 
are illustrated in Figure 10.2 for some particular values of a and b with n 
equal to 4, 6 and 8. 

In all cases the approaching waves become unbounded as / — > — | in 
region II, and as ^ — > — | in region III. 
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1 1 1 1 1 

1 H 1 1 1 1 1 

1 



(c) 6^-2, a= 2 -^2(1-^). (d) 6 = - 2, a = 2 + ^2(1 - 

Figure 10.2 Scaled wave profiles of approaching waves in the Ferrari-Ibanez- 
Griffiths solution for various values of a and b, and with n taking the values 4, 
6 and 8. 

It can also be shown that this family of solutions has the same sin- 
gularity structure as the class of Szekeres solutions as described in Sec- 
tion 9.3 and to which it reduces when b = 0. In region IV, there is a scalar 
polynomial curvature singularity on the space-like surface t = 1 on which 
i/j — 7r/2 and f + g — 0. It may be noted that when 2 < n < 4, the initial 
boundaries u = and v = of region IV contain a distribution-valued 
singularity. These initial boundaries are regular when n > 4. As in all 
other solutions, it may also be noted that there are non-scalar curvature 
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singularities in the initial regions II and III on the hypersurfaces on which 
/ = — I and 9 = respectively. 

10.4 The soliton solution of Ferrari and Ibanez 

Ferrari and Ibanez (19876) have also used a generating that is familiar in 
the study of solitons to obtain a solution of (10.2) in which 

V = -2kQoit)Po{z) - |alog(l - 1'^) - |alog(l - z'^) (10.37) 

where k = 1 and a is an arbitrary constant.^ This was initially obtained 
using the inverse scattering technique in which the initial 'seed' solution 
is taken to be the Kasner or Stoyanov solution given by (10.20) or (10.29). 

This solution can immediately be seen to be a generalization of the 
Khan-Penrose solution (10.18) by the addition of the solution (10.20). It 
necessarily contains impulsive wave components. 

It was pointed out by GriflBths (1987), however, that this solution 
may easily be generalized by treating the parameter k in (10.36) as 
another arbitrary constant. In this way, this family of solutions may 
also be considered to be a generalization of the Szekeres solutions with 
k\ = k2 = k. As in the previous example, this generalization then permits 
the approaching waves to have a continuous wave front for appropriate 
values of k and a. It is this more general class of solutions that is consid- 
ered in this section. 

It is possible immediately to integrate equations (10.23) to obtain 

S = const + k'^ log(t2 - z^) - la^ log(l - z^) 

- |(a - 2kf log(l - t) - l(a + 2kf log(l + t). ^^^'^^^ 

This contains the term k^ log(t^ — z'^) which, as in (10.33), is required to 
ensure that is continuous across the boundaries of region IV. If the 
leading terms in the power series for / and g take the form 

/ = I - (ci«)" + . . . , ^ = I - ic2vr + ■■■, n>2, (10.39) 

then the required boundary conditions are satisfied if 

P=2(l-l/n). (10.40) 



As well as putting A; = 1, Ferrari and Ibanez find it convenient to use the 
parameters Si = |(1 - a) and S2 = + a). 
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In this case it is remaxkable that the junction conditions place no con- 
straint on the parameter a. This feature is, in fact, related to a general 
result that will be described in Section 12.1. 

When evaluating expressions for the Weyl tensor, it is convenient to 
use the coordinates ip and A defined by (10.10) (it may be noticed that 
these differ from those used by Ferrari and Ibanez). The scale-invariant 
components of the Weyl tensor are given by 

_ g'^ f 2k{l - P) cos^ X 3k{ak -smip)cos'^ X 

^0 — ~ „^„2 „/ 2 \ ( _• _3/.;. 



cos^ ij; cos^ A \ sin^ (i/j - A) sin^ {ip - A) 

3k(l — a^)cos A a , n, 
+ sin(^-A) -z'l-"') 

,po ^ fs' ( k{k-asm ip) cos^ A _ l^_^2x\ 
^ cos^ Vcos^ A \sin(V' — A) sin(7/; + A) 4 J 
^o__ f'^ /2k{l-k^)cos^X 3k{ak-smi;)cos^X 
^ cos2 i/j cos2 A V sin^ {ip + A) sin^ {ip + A) 

3/c(l — a^) cos A a. 2^ 
^ sin(V; + A) 2^ ~ ' 

(10.41) 

For the sake of later discussion, it is found to be appropriate to choose 
the scale functions A and B defined by (6.2) and (6.14) to be 



(10.42) 



B = 



f sin'='/2(V.-A)'' 

(cos COS A) sin'^' {i; + A) 

7'(1 + sinV')'=(^+")/2(l -sinV')'=('=-")/2 

(/ + ^)V^sin^V^(V;-A) 
9' sin'='/2(V' + A) ^ 

(cos v' COS A)(^ -"'V4 _ A) 

- ~ g'(l + ^i-a_'^)k{k+a)/2(^l _ ^)fe(fe-a)/2 " 

The negative signs are due to the signs of /' and g' . 

From (10.41) it can be seen that, in general, this solution has the 
same singularity structure as the Szekeres solution. There is a curvature 

singularity on the hypersurface given by/ + (7 = 0or'i/' = 7r/2. 

There arc, however, exceptional cases which occur when a = ±1 and^ 



^ Since V can always be replaced by - V, there is no loss of generality by 
assuming that k is always positive. 



10.5 The degenerate Ferrari-Ibanez solutions 



81 



k = 1. In these cases 



3 



(10.43) 



(1 ± simp)^ 



which satisfies the condition 



*0*4 = 9*2 



(10.44) 



which imphes that the space-time is of algebraic type D (see Kramer et 
al. 1980, or Chandrasekhar and Xanthopoulos 19866). These particular 
degenerate cases will be analysed in more detail in the next section. 

When n — 2 and k = 1, the approaching waves contain an impulsive 
component. For 2<n<4, l<A;^<3/2, the null boundaries of 
region IV contain distribution-valued singularities. These boundaries are 
regular when n > 4, 3/2 < A;^ < 2. The approaching waves have a step 
wavefront if n = 4, and the wavefront is continuous if n > 6. 

Also, as in the Szckcrcs solutions, there are non-scalar curvature sin- 
gularities in regions II and III when / = —1/2 and g = —1/2 respectively, 
for all values of a and all permissible values of k. 

A generalization of the class of solutions described in this section has 
been obtained by Tsoubelis and Wang (1989). This has been obtained by 
putting 



V = -ialog(l - t2)(i - z^) - hQo{t)Po{z) - b2Po{t)Qo{z). (10.45) 



By a comparison with (10.37) and (10.19), this can be seen to be a general- 
ization both of the above Ferrari-Ibafiez solution and also of the Szekeres 
class of solutions described in Chapter 9. Its properties can reasonably 
be inferred from those of these two subclasses, as has been confirmed by 
Tsoubelis and Wang. 

10.5 The degenerate Ferrari— Ibanez solutions 

Consider now the special cases of the above Ferrari-Ibanez solution in 
which k = 1 and a = ±1, in which the Weyl tensor components are given 
by (10.43). The properties of these cases have been further described by 
Ferrari and Ibaiiez (1988). 

In the case when a = —1, k — 1, there is again a curvature singularity 
when i/j = Tv/2. However, when a — 1 and k = 1, the space-time appears 
to be regular for all < ip < 7r/2, and the singularity caused by the 
mutual focusing of the two waves appears to have been removed. It is, 
therefore, appropriate to consider this particular case in more detail. 
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Since A; = 1 in these cases, it is possible to scale the null coordinates 
such that 

f = l-u'^e{u), g = l-v^Q{v). (10.46) 
When a = 1, the metric functions in region IV are then given by 

1 /I — u^/l — — v^/l — 

e 



e 



1 - u'^ - v'^ yi + uVl^^ + v^/T^^ ) (10.47) 



Vi — tt^Vi — 

and the Weyl tensor components are then given by 

*o = ^ , S(v) ^ o G(^;) 

*2 = — ^ e{u)&{v) (10.48) 

(1 + ttVl -v^ + vVl - u^) 

1 3 

*4 = , 5(u) - o Q(u). 

As expected, these indicate that the approaching waves are impulses fol- 
lowed by continuous components. Of greater significance, however, is 
the fact that in this case there are two point singularities at the points 
u = 0,v = 1 and u = l,v = 0. The existence of the singularities at 
these points implies that the lines u = 1 and v = 1 in regions II and III 
act as 'fold singularities' like those described in Section 8.2. These lines 
therefore form boundaries to the space-time in these regions. 

It is convenient now to return to the time-like and space-like coordi- 
nates V' and A. In the degenerate case when a = 1 the line element in the 
interaction region takes the form 

d.^ = (1±^ (dV,^ - dA^) - ( l^^) _ eos^ A(l + sin 
2 \1 + smip J 

(10.49) 

To analyse this particular case further, consider the change of vari- 
ables 

r = 1 + sinih, r = V2x, 

(10.50) 

^ = 7r/2 - A, = V2y. 
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With this, the hne element (10.49) takes the form 

2ds2 = 1^1 - dr^ - -jy^dr'^ - ^^(d^^ + sm^ (10-51) 
which may immediately be recognized as the Schwarzschild metric with 

171 = 1. 

With < '0 < 7r/2, however, we have 1 < r < 2 which is normally 
understood as the region inside the horizon. In addition, since •0 increases 
from zero in the interaction region, this region must correspond to the part 
of the Schwarzschild space-time indicated in Figure 10.3(b) which is inside 
the initial horizon. However, it may also be noticed that the coordinate (p 
which represents the axial coordinate in the Schwarzschild solution covers 
the entire range —oo < < oo in this case. 




(a) X = const, y = const. (b) A = const, y = const. 

Figure 10.3 The interaction region of the degenerate solution viewed in two 
different coordinate planes. The plane (a) is the u,v or i/jjX plane as in previous 
figures. Diagram (b) represents the i/jjX plane and is equivalent to the familiar 
Kruskal-Szekeres diagram for the Schwarzschild space-time. 



Prom Figure 10.3, it appears that the approaching waves collide at 
the surface '0 = and reach a horizon at V' = 7r/2. There is, however, 
no a priori reason why the coordinate il^ should not be continued beyond 
7r/2. If such a continuation is possible, this would indicate that the waves 
would continue through the horizon until they finally end in a curvature 
singularity at i/j = 37r/2. 
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The other degenerate case in which a = — 1 can also be transformed 
to a Schwarzschild space-time but, in this case, the interaction region cor- 
responds to the upper region inside the horizon and all trajectories end in 
the future singularity. This case, together with the above transformation, 
has also been described by Yurtsever (1988a). 

It is instructive also to express the metric in Kruskal-Szekeres form. 
This can be achieved from (10.49) using the transformation 



u = 



-VI - sin ^ e(i+''" '^^/^ e^/^ for < V < 7r/2 

-I- VI - sin^A e(i+^'" '^^/^ e^/^ for n/2<i;< 3tt/2 

-VI - sin 6^^+^^" e"^/^ for < V < 7r/2 

+V1 -sinV' 6^1+^'"'^)/^ e-^/^ for 7r/2 < V < 37r/2 



(10.52) 



where u and v are null coordinates in the ip-x plane. With this, the line 
element becomes 

2ds2 = e-^^+^^^^)/2j~j~ _ Q ^ sinV')^(d^2 + sin^ Odcf)^) 

(1 -I- smip) 

(10.53) 

where i/j is given by 

{l-smi;)e^/^^^ = uv. (10.54) 

In these coordinates the structure of the solution in the vicinity of the 
horizon is clearly shown to be regular. 

It is also convenient to consider this solution in terms of alternative 
null coordinates^ u' and v' such that, in the interaction region 

J 1 ■ 2 / 1 -27/ 

1 = 7:— sm au , <7 = o — sm ov 

, , , (10.55) 

ip = au +bv , X = au —bv 

where the new constants a and b are the strengths of the approaching 
gravitational shock waves. With these coordinates it can be seen that the 
Schwarzschild mass parameter, which is unity in the line elements (10.51) 
and (10.53), is here related to the amplitudes of the approaching waves 

by 

1 

m=^=. 10.56) 
V2^ 



^ These coordinates were initially found useful by Bell and Szekeres (1974) 
and will be introduced here in Section 15.1. They have been used to analyse 
the degenerate Ferrari-Ibanez solution by Hayward (1989a). 
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It follows from this that stronger approaching gravitational waves would 
produce a shorter proper time between the collision and the subsequent 
horizon. They would also give rise to a smaller Schwarzschild mass pa- 
rameter, and hence greater curvature on the horizon. 

When trying to relate the global structure of this solution to part of 
the Schwarzschild space-time, it should be emphasised that for this class 
of colliding plane waves the coordinate (f) is not periodic but covers the 
full range — oo < < oo. Topological singularities occur on the lines 
^ = and ^ = TT on the horizon ip — 7r/2 (or r = 2). 

This global structure of the solution has been determined by Hayward 
(1989 a). He has noted that the above solution also forms a covering space 
of part of the Schwarzschild white hole with a quasiregular covering space 
singularity along the polar axis around which the space-time is wound. He 
has also suggested an alternative extension of the Schwarzschild exterior 
inside the black hole. 

In addition, he has suggested some possible extensions for the non- 
unique region beyond the horizon including one that is time symmetric. 
Refering to Figure 10.3(b), the opposing waves in this case collide on 
the hypersurface ip = 0. The interaction region is then that part of the 
Schwarzschild space-time inside the initial horizon as indicated, which ex- 
tends to the horizon as the focusing hypersurface. An extension through 
this surface may be taken to be the familiar two parts of the exterior 
Schwarzschild space-time. This can be continued to the subsequent hori- 
zon. A further extension through this horizon is then possible up to the 
next surface on which t/j = or r = 1 at which the space-time again splits 
into two separating gravitational waves which are the time reverse of the 
initial approaching waves. 

At the conclusion of this section, it may be noted that Ferrari and 
Ibanez (1988) have shown that the shear-free principal null congruences 
associated with this type D space-time do not focus on the horizon. How- 
ever, these congruences are in different planes from the congruences on 
which the two wave components propagate, and for which the contraction 
and shear clearly become unbounded as •0 — > 7r/2. 

10.6 An odd order solution 

All of the solutions considered so far have involved only Legendre funtions 
of even order. A solution with odd order functions has been presented by 
Griffiths (1987)^. This has 



This solution with o = 2 has subsequently also been given by Li (1989). 




(10.57) 
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With this function, equations (10.23) can immediately be integrated to 
give 

S = const + log(t^ - z'^) - log(l - z^) 



As in (10.32-34), it can be clearly seen that this expression contains terms 
of the form (7.10). It may thus be concluded that, if the leading terms 
in the expansions for the functions / and g are given by (10.34), then 
the boundary conditions that are required for the solution to describe 
colliding plane waves are satisfied if 



As in previous examples, it can be shown that the approaching waves 
have an impulsive component if n = 2. Also the wavefront is unbounded 
if 2 < n < 4, has a step if n = 4, is continuous but not smooth if 
4 < n < 6, and is smooth if n > 6. The expressions for the Weyl ten- 
sor components are rather complicated, and there is nothing particularly 
remarkable about this solution. 

In the previous examples discussed in this chapter, the two waves 
approaching each other are identical. The component \l/o in region III 
has an identical form to ^'4 in region II but with the null coordinate 
V replacing u. Similarly, the expressions for V in the two regions can 
be related by the same replacement. In this case, however, since odd 
functions are being used, a change of sign is also included. Thus, V{u) 
and \E'4('u) in region II are replaced by —V{v) and — \l/o(^) in region III. 
The approaching waves still have colinear polarization but, in this case, 
their amplitudes are opposite rather than the same. 

10.7 The second Yurtsever and the Feinstein— Ibanez so- 
lutions 

For colliding plane wave solutions, we are considering space-times with 
two space- like Killing vectors. Such space-times have also been consid- 
ered in the context of cosmology, where certain vacuum inhomogeneous 
cosmologies satisfy the same field equations as the interaction region for 
colliding colinear gravitational waves (Gowdy 1971), although different 
boundary conditions are appropriate. In fact the vacuum Gowdy cos- 
mologies can be considered to represent closed universes built from op- 
posing plane gravitational waves. The similarity between these solutions 




(10.58) 



o2/8 = 1 - 1/n. 



(10.59) 
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has been noted by Feinstein and Ibanez (1989), who have used known 
cosmological solutions to obtain a new class of colliding plane wave solu- 
tions. In this section, the Feinstein-Ibahez solutions are presented. The 
electromagnetic Gowdy cosmologies will be described in Section 17.3. 

It is appropriate at this point to consider the alternative system of 
coordinates defined by 

i=/ + 9 = VT^Vw ^^^^^^ 

z = f - g = -tz. 

It may be noted that the new coordinate t is a decreasing, or past pointing, 
time-like coordinate, and that the singularity in region IV occurs when 
i = 0. These coordinates have also been used by Yurtsever (1988c)^. 

In this case, the line element (10.1) for colliding waves with aligned 
constant polarization can be rewritten in the diagonal form 

ds^ = ^ (dp - dz2) - t (e^dx2 + e-^d2/2) (10.61) 
and the main vacuum field equation (9.3), (10.2) or (10.13) becomes 

V+lv-V"^0 (10.62) 
and the subsidiary equations (10.3) become 

S = \t (y^ + V'^^ , S' = -iVV (10.63) 

where the dot and prime denote derivatives with respect to i and z re- 
spectively. 

The above equations (10.62,63) are exactly those for the vacuum 
Gowdy cosmologies that have been considered by many authors.^ The 
general solution of (10.62) can be expressed as a line integral. However, 
for the class of Gowdy cosmologies, it is found to be convenient to consider 
the class of solutions given by 

V = - alogi + Li{A^ cos[u;{z+au,)]Jo{u;i)} 

+ L2{B^ cos[uj{z + P^)]YQ{ujt)} + Y,di cosh-i (^7^) ^^^'^^^ 



^ In the notation used by Yurtsever t and z are replaced by a and /3 
respectively. 

^ See for example Centrella and Matzner (1979), Moncrief (1981), Carmeli, 
Charach and Malin (1981), Carr and Verdaguer (1983), Ibanez and Verda- 
guer (1983), Adams et al. (1982), Feinstein and Charach (1986), and Feinstein 
(1987). 
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where Li{ } and -^2{ } are arbitrary linear combinations of the terms in 
curly brackets, including the Fourier-Bessel integrals of the form 



Jo(a;f) and lo(<^^) are Bessel functions of the first and second kinds of 
zero order, and A^^, B^^, a^o and (3^^ arc sets of arbitrary constants. 

The first term in (10.64) is clearly identical to the solution (10.20). 
It is the term that, on its own, gives rise to the Kasner solutions. The 
set of terms included in Li have regular behaviour as t — > 0. It may 
be noted that this property is also shared by the combination of terms 
'^anPn{t)Pn{z) that are contained in (10.16). The set of terms included 
in L2, however, are badly behaved as t ^ 0. These are the terms that 
are considered to induce chaotic behaviour in these cosmological models. 
The last term in (10.64) is identical to (10.31). This contains the so-called 
gravitational soliton components. 

Yurtsever (1988 c) has considered the asymptotic structure of those 
solutions which only contain the combinations Li{ } and }, i.e. in 
which a — di — 0. He has shown that these solutions are asymptotic to the 
inhomogeneous Kasner solutions as the singularity f = is approached. 
He has also given explicit expressions which relate the asymptotic Kasner 
exponents along the singularity to the initial data specified along the wave 
fronts of the incoming colliding plane waves. It follows from this analysis 
that the focusing hypersurface t = is a curvature singularity except in 
the special case in which one of the Kasner exponents is zero. This special 
case of the degenerate Kasner solution is flat and, in this case, the focus- 
ing hypersurface is a Killing-Cauchy horizon across which space-time can 
be extended. It is reasonable to conclude from this that, although there 
are colliding plane wave space-times which contain a Killing-Cauchy hori- 
zon rather than a space-like curvature singularity, these space-times are 
unstable against small perturbations of the initial data and that 'generic' 
initial data always produce space-like curvature singularities. 

Feinstein and Ibahez (1989) have considered the family of solutions 
for which the combination L2 is zero. They have shown that, for this case, 
a curvature singularity does not develop in region IV as t — > provided 
that 



They have also shown that, in this case, the necessary boundary condi- 
tions (7.15) are satisfied provided there are at least two solitonic terms 




(10.65) 




(10.66) 
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with constants satisfying 

ci = -1, c^i^ = 2(1 - i/m) 

\ ^ ^ ^ (10.67) 
C2 = l, = 2(1 - l/na) 

using the notation of (7.11). (The unfortunate repeated use of constants 
ci and C2 in this equation should not cause confusion.) It may thus be 
noted that these two soHtonic terms provide the conditions for continuity 
on the two different null boundaries of region IV. For continuity it is also 
required that ni > 2 and n2 > 2, so di and ^2 are constrained to the 
range 

1 < |c^i|, |(i2| < V2. (10.68) 

It follows from this and (10.66) that a 7^ 0, and so the first term in (10.64) 
must necessarily be included in these solutions. 

It may also be noted that the solutions given in the notation of 
previous sections by 

V = - |alog(l - t^) - |alog(l - z^) 

n i \ V V / 

where ci and C2 are given by (10.67), similarly do not contain curvature 
singularities on the hypersurface f + g = provided the constants di and 
a are constrained by (10.66) and (10.68). 

It is also of interest to note that the two solitonic terms that provide 
the continuity across the boundaries of the interaction region are in fact 
identical to the two separate terms in the Szekeres solution (9.4). This 
may be observed by noting that 

cosh""*^ I , I = cosh~^ 

= 2tanh~ 




-log I ^ 7=1 (10.70) 

= -Qo{t)Po{z) - Po{t)Qo{z) 
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and similarly 

cosh-i (^_J==i|_^ = -Qoit)Po{z) + Poit)Qoiz). (10.71) 

It follows from this that even the Szekeres solution, which is every- 
where at least C^, can be adapted by the inclusion of a suitable multiple 
of the term (10.20) to provide a solution without a curvature singularity in 
the interaction region. Thus, the occurrence of a quasiregular singularity 
in region IV can have nothing to do with the relaxation of the continuity 
conditions across the boundaries. 

The class of solutions included here contain an arbitrary number of 
parameters. It can be shown that they are algebraically general in the 
interaction region. They contain the usual coordinate singularity on the 
hypersurface f + g — but, for this class, this is not a curvature singu- 
larity. Feinstein and Ibanez have shown that the solution is extendable 
across this surface. However, the extension is not unique. 



10.8 The first Yurtsever solutions 

For colliding plane waves with aligned linear polarization, it has been seen 
that the main field equation (10.2) is linear. Various classes of solutions 
have already been obtained by separating the variables of this equation 
in a number of different ways. A further class of solutions in which the 
main equation is separated in yet another way has been given by Yurt- 
sever (1988 a). These solutions were originally obtained indirectly by first 
considering the Weyl solutions for stationary axisymmetric space-times. 
However, they will be presented here in a more general way that is also 
more appropriate in considering colliding wave solutions. 

It is appropriate here to start with the main equation in the form 
(10.62) using the variables i and z defined by (10.60), and then trans- 
forming it by putting 

f=i/sinh?7, z = i'COsh.r] (10.72) 

where the parameters v and rj are not necessarily real. Indeed, real values 
of these variables only cover part of the interaction region as will be 
clarified later. With these parameters, the main equation (10.62) becomes 

+ -K - 4 cothrjVr^ - J-F = (10.73) 
and this has a series of separable solutions of the form 

^ = ^ (cnZ^" + dnV-""-^) {anPnicOshrj) + bnQn{cOshT])) (10.74) 
n 
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f+g=o 




Figure 10.4 In the interaction region IV for the Yurtsever solutions, the vari- 
ables u and cosh rj are imaginary in the initial region marked A that immediately 
follows the collision. This region is followed by the regions marked B in which 
these variables are real. 



where Pn{cosh.r]) and ^^(coshry) are Legendre functions of the first and 
second kinds. 

We must now reconsider the parameters of these solutions. As defined 
by (10.72), u and coshry are imaginary when the product fg is positive. 
Since both / and g decrease from | towards — | in the interaction region, 
it is appropriate to divide this region into the subregions A and B in 
which the product fg is positive and negative respectively as indicated in 
Figure 10.4. Explicitly we may put 



In A : v = 2i\fjg., cosh 77 = 
In B : V = 2-^— /^f, coshry = 



ij^ (10.T5) 



It is clear that u is zero midway along the two boundaries between 
region IV and regions II and III. The inverse functions of v that appear in 
(10.74) must therefore be excluded for colliding plane wave solutions since 
they must be regular along these boundaries. It is therefore appropriate 
to express general solutions in the form 

v-.i:.w-(.„.„(|g))..„.„(|g))) 

(10.76) 



92 



Other vacuum solutions with aligned polarization 



This solution will from now on simply be quoted in the second form since 
either equation formally includes the other. 

It may be noted that these solutions are continuous on the bound- 
aries between the subregions A and B. Also, the terms involving the Leg- 
endre functions of the first kind are regular on the focusing hypersurface 
/ + ^ = 0, while those involving the Legendre functions of the second kind 
are singular on this surface. 

To be precise, in the original class of solutions given by Yurtsever 
(1988a), V is taken in the form 



which can be seen to be generalizations of the degenerate Ferrari-Ibanez 
solution described in Section 10.5. They are therefore distortions of the 
Schwarzschild black hole solution in the interaction region, and corre- 
spond to interior Weyl solutions for static axisymmetric space-times. In 
terms of colliding plane waves, these solutions all involve approaching 
waves with initial impulsive components. Some particular examples have 
been described explicitly by Yurtsever (1988 a). 

Clearly these solutions can easily be generalised by modifying the 
coefficients of the first terms of (10.77). By choosing these appropriately, 
it is possible to remove the impulsive components from the approaching 
waves. 

10.9 Further explicit solutions 

To obtain solutions describing the collision and interaction of plane waves 
with aligned linear polarization it is necessary first to solve the main field 
equation, which may be expressed in terms of different coordinates in any 
of the forms (9.3), (10.2), (10.13), (10.62) or (10.73). 

Any number of further explicit solutions may easily be generated 
using the methods described in the previous sections of this chapter with 
different combinations of particular solutions of the main equation (9.3), 
(10.2), (10.13), (10.62) or (10.73). It is appropriate here simply to list the 
various possibilities. 

First there are the solutions (10.16) which have been obtained by 
separating the variables in (10.13). These involve Legendre functions of 
the first and second kinds 



Vl^J2 [(^nPnit) + qnQn{t)]Pn{z) + [PnPn{t) + bnQn{t)]Qn{z) 
n n 

(10.78) 



V = -2Qo{t)Po{z) + 



|log(l-t2) + ilog(l-^2) 




(10.77) 
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where a„, g„, Pn and bn are series of arbitrary constants. It may be noted, 
of course, that non-integer values of n may also be included. 

There are also the particular solutions (10.20), (10.8) and (10.22) 
that can be re-expressed in the form 

V2 = -alog(/ + 9) + Y. , 

.E*--.(^) 

where a, A^, ai, ci and di are arbitrary constants. 

Then there are the solutions of (10.62) that are contained in (10.64) 
and have previously been included in Gowdy cosmologies, namely 

^3 = ^A^cos[u{z + a^)]Jo{ut) + ^B^cos[u{z + l3^)]YQ{ui) (10.80) 

where Jo{uji) and YQ{ui) are Bessel functions of the first and second kinds 
of zero order, i=f + g, z = f — g, and A^J, B^, and P,^ are sets of 
arbitrary constants. 

It may be noted that the solutions (10.80) have been obtained by 
separating the variables in (10.62), taking only solutions that are periodic 
in z. This constraint is not required for colliding plane waves and an 
additional set of aperiodic solutions is given by 

14 = ^ Ai_j cosh[ci;(5 + a^)]Io{uji) + ^ B^j cosh[a;(5 + (3aj)]Ko{uji) 

(10.81) 

where Io{u3t) and ^0(0;^) are modified Bessel functions of the first and 
second kinds of zero order. 

Finally, there are the separable solutions of (10.73) which are non- 
singular on the initial boundaries, namely 




(10.82) 



which again involve Legendre functions of the first and second kind. 

When considered as infinite series, these different forms may simply 
be regarded as different representations of the same class of solutions. 
However, when looking for particular explicit solutions it is necessary to 
consider only a few terms. In this way, further exact solutions for colliding 
plane waves may be obtained by combining particular components from 



(10.79) 
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any of the series Vi, V2, V3, V4 and and by choosing the arbitrary 
constants such that 

It,, [C2 - 9)Vn = 2/C2, lim [(I - f)V^] = 2h (10.83) 

where ki and /c2 are constrained to the range (7.13). If this condition is 
satisfied, then it remains only to integrate the subsidiary equations in the 
form (10.3) or (10.63) to obtain the remaining metric function S or M. 



